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Abstract We report on an experimental investigation of the tumbling of mi- 
crorods in the shear flow of a microchannel (dimensions: 40 mm x 2.5 mm x 
0.4 mm). The rods are 20 |am to 30 jam long and their diameters are of the 
order of 1 \im. Images of the centre-of-mass motion and the orientational dy- 
namics of the rods are recorded using a microscope equipped with a CCD 
camera. A motorised microscope stage is used to track individual rods as 
they move along the channel. Automated image analysis determines the po- 
sition and orientation of a tracked rods in each video frame. We find different 
behaviours, depending on the particle shape, its initial position, and orien- 
tation. First, we observe periodic as well as aperiodic tumbling. Second, the 
data show that different tumbling trajectories exhibit different sensitivities to 
external perturbations. These observations can be explained by slight asym- 
metries of the rods. Third we observe that after some time, initially periodic 
trajectories lose their phase. We attribute this to drift of the centre of mass 
of the rod from one to another stream line of the channel fiow. 

Keywords Microrods • microchannel • Jeffery orbits • tumbling 



1 Introduction 



It was shown by Jeffery (1922) that small axisymmetric rods in a viscous 
shear flow align for the most time with the flow direction, but that their 
symmetry axes periodically and rapidly turn by 180 degrees. This motion is 
referred to as 'tumbling' in the literature. Experimental studies of tumbling of 
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particles in flows have been performed for a long time, see for example ( Gold- 



smith and Mason, 1962). More recently, Kaya and Koser (2009) performed an 



experiment analysing the tumbling of E. coli cells in microchannel shear flows. 
They characterised the tumbling by fltting periodic orientational trajectories 
to short experimental time series. A similar approach was adopted by |Mishra| 
[et al.| ( ,2012, ) to describe the tumbling of microrods in a micr ochannel flow . 
The periodic solutions referred to above were flrst obtained by Jeffery (1922) 
and are commonly referred to as 'Jeffery orbits'. 



Hinch and Leal ( 1979| ) have shown by theoretical analysis that slightly 



asymmetric rods also tumble, but that their orientational motion is in general 
not strictly periodic. Under certain circumstances t he authors p redict that 
the orientational dynamics may be 'doubly periodic'. Yarin et aL] (p!997 ) refer 
to this motion as 'quasi periodic'. Their numerical results show that the 
tumbling can be chaotic (Ott, 1993). In order to experimentally distinguish 



periodic from aperiodic tumbling due to asymmetry, it is necessary to observe 
long sequences of flips. A confounding factor is rotational and centre-of-mass 
diffusion: the centre of mass of an advected rod may diffuse to neighbouring 
stream lines of the flow. This causes aperiodicity. But noise can also directly 
affect the rotational degrees of freedom, leading to random tumbling. 

In order to disentangle these effects it is necessary to follow the orien- 
tational dynamics for many flips. The rods must therefore be tracked for 
long distances. To achieve this we have improved an existing experimental 
setup (Mishra et al. 2012) in two ways. First, we have automated the image 
analysis of the empirical data. Second, with the new setup it is possible to 
periodically revert the channel flow. This allows us to record longer tumbling 
sequences. In an ideal experiment the orientational dynamics is expected 
to retrace its trajectory upon reversal of the flow. This is a consequence of 
the time-reversal invariance of the Stokes equation governing low Reynolds- 
number flow. Our setup thus enables us to quantify the sensitivity of the 
observed tumbling motion to perturbations. In the following we report and 
discuss experimentally observed tumbling trajectories. First, we see periodic 
and aperiodic tumbling. Second, different tumbling trajectories exhibit differ- 
ent sensitivities to external perturbations. We argue that these observations 
can be explained by slight asymmetries of the rods. Third, at very long times 
initially periodic trajectories lose their phase. We attribute this to drift of 
the centre of mass of the rod from one to another stream line of the channel 
flow. 

We conclude this introduction by briefly commenting on the relevance of 
the questions addressed in this paper. Jeffery's periodic tumbling solutions 
in shear flows form, together with orientational diffusion, the basis of many 
st udies of the rheology of suspensions o f non-spherical particles, as outlined 
by Hinch and Leal (1972). Petrie (1999) has given an overview over the rhe- 



ology of flbre suspensions and the use of Jeffery's and diffusion theory in this 
fleld. Furthermore, pattern formation by non-s pherical particles in ran dom 
flows was inves tigated using J effery's theory by Wilkinson et al. ( 2009 ) and 
Bezuglyy et al.| ([2010) (see also [Wilkinson et_aL (2011)), identifying singular- 



ities in the orientational patterns of rheoscopic suspensions, and explaining 
how Jeffery's periodic solutions determine rheoscopic visualisations of flows. 
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Last but not least, the tumbling and alignment of small rods in turbulent 
flows has recently been intensively investigated, both experimentally, by sim- 
ulations, and theoretically. We refer to the art icles i n this special issue of Acta 
Mechanica, as well as to (Parsa et al. , 2012) and ( [Wi lkinson and Kennard 



[2012) and articles cited in these two papers. Simulations and the theoreti- 
cal treatments of the tumbling of small rods in turbulent flows are based on 
Jeffery's equation of motion. 



2 Materials & methods 



In summary, the microrods are produced by means of a liquid-liquid disper- 



sion technique (Alargova et al. , 2004). They are introduced into a micro- 



channel system produced in PDMS. The flow is pressure-driven, using a sy- 
ringe pump. The flow direction is periodically reverted. The fluid is a 2:1 
mixture of glycerol to water, so the viscosity is high, and the Reynolds num- 
ber is low. The centre-of-mass motion and orientation of the microrods are 
observed using a 20 X microscope equipped with a CCD camera. A motorised 
stage allows to track individual rods over long time periods. The resulting se- 
quence of frames is analysed by an automated image-analysis algorithm pro- 
grammed in MATLAB, resulting in the orient ational trajectories displayed 
in Figs. [5] to [T) 



2.1 Experimental methods 



Particle synthesis. As in our earlier experiments ([ Mishra et al.[ [2012 ), the 
polymer microrods are prepared by a liq uid-liquid dispersion technique us- 
ing the protocol of Alargova et al. (2004). This method produces rods with 
lengths of 10 \im to 100 |xm, with typical aspect ratios of the order of 10:1. 
The resolution of the microscope does not allow us to determine whether 
the rods are symmetric or not, that is, whether they have perfectly circular 
or slightly elliptical cross sections. Inspection under an optical microscope 
shows that the shorter rods are straight. Some of the longer rods are curved. 
In the experiments described below, rods in the shorter range of the span 
given above are used. 

Design and fabrication of microfluidic channels. The microfluidic device 
used in this experiment is shown schematically in Fig.jl] The channel (40 mm 
long, 2.5 mm wide, and 0.4 mm d eep) is pro duced in PDMS. The procedure 
is described in detail by Mishra et al. (2012). 

Optical system and tracking. A sketch of the experimental setup is shown 
in Fig.[T] The optical setup is based on a Nikon Eclipse Inverted Microscope. 
The channel is placed on a motorised stage, making it possible to track the 
rods as they follow the flow in the channel. The rods are observed using 
a 20X microscope objective (NA=0:28). Images are recorded with a CCD 
camera (Leica DEC 350EX, Switzerland) at a frame rate of 100 frames per 
second. As shown in Eig. [l] the images are observed through the cover glass. 
Since the magnification of the microscope is only 20X, the thickness of the 
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Fig. 1: Experimental setup (schematic) with coordinate system. The mi- 
crochannel is placed upon a motorised stage which allows the camera to 
study any part of the channel. The calculated flow rates in a cross-section 
of the channel are shown in Fig. [2] 



glass does not cause focal problems. The channel is illuminated through the 
transparent PDMS. 

Fig. [l] also shows the Cartesian coordinate system adopted in this paper. 
The X-axis is taken to lie along the flow in the channel, the 7/-axis lies along the 
optical axis of the microscope (the channel is 400 \im deep in this direction) . 
The z-axis, finally, lies along the width of the channel (it is 2.5mm wide). 
Thus the camera plane corresponds to the x-z-plane, and the channel cross 
section (Fig. [2]) corresponds to the y-z-plane. 

Channel flow. The flow in the microchannel is characterised by a very 
small Reynolds number. The rods are suspended in a 2:1 mixture of glycerol 
and w ater, corresponding to a kinematic viscosity ofu = 2-10~^ m^/s (Cheng 



2008). The flow speed is of the order of 100|am/s. Based on the smallest 
channel dimension (400 \im) this yields a Reynolds number of the order of 
Re = 10~^. The Navier-Stokes equations for the incompressible flow in the 
channel thus reduce to the viscous Stok es equation which c an be simply solved 
by a basis expansion. Closely following |Brody et al. (1996) we have computed 
the flow profile in a cross section far from the in- and outlets (so that we can 
neglect the x coordinate). We assume no-slip boundary conditions at the 
channel walls and a constant pressure gradient over the channel length. The 
resulting profile for our channel geometry is shown in Fig. [2] The resulting 
flow velocity scales as Ux ~ Ap/rj as a function of pressure difference Ap and 
dynamic viscosity r]. But the form of the profile is not affected by the values 
of Ap or T]. In Fig. [2] we observe that close to the centre of the channel, the 
flow gradient is oriented along the y-axis. In the following we refer to this 
axis as the 'shear direction'. The x-axis is termed 'flow direction', and the 
2;-axis 'vorticity direction'. 
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Fig. 2: Computed flow profile for the channel geometry used in experi- 
ments. Width {z) of channel is 2.5 mm and depth (y) is 400 \im (refer to 
Fig. [l] for coordinate system) . The intensity plot shows a cross-section of 
the channel, intensity indicating flow velocity in arbitrary units. Top and 
right panels show the corresponding flow- velocity profiles at y = and at 
z = 0, respectively. 

The rods analysed in this paper are between 20 \im to 30 \im long. The 
viscous time scale L"^ /u is of the order of 10~^s, which is much smaller than 
both the tumbling and the flow-reversal time scales. The disturbance flow 
around the inclusion is therefore expected to obey the quasi-steady Stokes 
equation, in spite of the motion of the rod. 



2.2 Data analysis 

As explained below, the rods often tumble aperiodically. To classify different 
dynamical behaviours requires long time series. Since one traversal of the 
channel takes approximately 5 minutes, large amounts of data must be anal- 
ysed. We have therefore designed a data analysis software (in MATLAB) to 
track the centre-of-mass and orientational motion of the suspended rods in 
an automated fashion. 

The raw data from the experiment consist of time-stamped positions of 
the motorised stage, together with frame sequences from the microscope. The 
analysis uses this information to estimate time series of the centre-of-mass 
position of a given rod, and the orientation of the unit vector n along its 
major axis: 

n = {nx.fiy.riz) . (1) 

The coordinate system is shown in Fig. [l] 

In the remainder of this section we briefly describe how the empirical data 
were processed. First, we summarise our algorithm for the image analysis. 
Second we explain how time is rescaled to account for variations in the flow 
velocity due to reversals of the direction of the flow and due to fluctuations 
in the pump pressure. 
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(a) 

Fig. 3: Three steps of the data analysis, (a) the particle is located within 
the video fram e, (b) edges are detected with an edge-detection algorithm 
(Canny, 1986| ), (c) an ellipse is fitted to the resulting edge pixels. The ori- 
entation of the major axis of the ellipse is used to estimate orientation of 
the rod in the camera plane. 



2.2.1 Image analysis 

The automated image- analysis software tracks the centre of mass and the 
orientation of a given rod as it moves in the channel. We separately discuss 
two steps: detecting the position and orientation of the rod, and estimating 
its length. 

Finding the rod position and orientation. To begin with, static image noise 
is removed as follows. For each pixel, the time-averaged intensity is computed, 
and for each time, the average intensity of the frame is calculated. Then, for 
each pixel, the time-averaged intensity is subtracted from the instantaneous 
intensity, and the instantaneous frame average is added. 

To locate the rod in a given frame, the software requires an initial guess 
of the centre-of-mass position. The frame rate is sufficiently high so that 
the tracked rods move only slightly between subsequent frames. The initial 
guess is therefore taken to be the last-known rod position (for the first frame, 
manual user input is required). 

Once the location of the rod within a given frame has been estimated, a 
smaller window is cropped around the approximate position of the particle. 
An example is shown in Fig. |3a| A standard implementation of the edge- 
detection algorithm suggested by Canny] C1986) is used to find the points 
on the edge of a rod sufficiently contrasted from its surroundings, as shown 
in Fig. |3b| Once the positions of the edge pixels of the r od in question are 
known, a least-squares fit of an ellipse to its edge pixels ( Halif and Flussei] 
1998) is used to find the shape and orientation of the projection of the rod 
in the camera plane (Fig. 3c). The n^-component is not directly observable, 
but can be computed from rix and if the length of the rod is known. 

Estimating the rod length. It is impossible to tell the length of the rod from 
the length of its projected image in a single frame. However, the rod spends a 
significant amount of time aligned with the flow in the x-2;-plane, where the 
projected length corresponds to the rod length. We therefore recorded time 
series of projected rod lengths and computed the distribution of projected rod 
lengths (Fig.|4|. The position of the maximum of the distribution determines 
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Fig. 4: Histogram of observed projected rod lengths in the x-z-plane fol- 
lowing a single particle through the flow. The rod spends a long time 
nearly aligned with the flow direction. This results in the peak of the his- 
togram. The tail to the left corresponds to observations while the rod is 
turning around, producing a shorter projection in the camera plane. The 
peak is broadened by measurement noise due to the finite resolution of the 
camera. The rod length is taken to be the position of the peak. Data corre- 
sponds to rod #1 (see Section [s]). 

the rod length. We note, however, that the empirical histogram is subject to 
noise (as a result of the finite resolution of the microscope). This may cause 
a bias in our estimate of the length. 

2.2.2 Time rescaling 

As explained above, we revert the flow direction in order to test the sensitivity 
of observed orientational motion to perturbations. We have found that small 
perturbations may have significant effects. It is therefore important to revert 
the flow as smoothly as possible. This means that the flow velocity changes 
substantially over a significant amount of time. Another source of changes 
in flow speed Ux is noise in the form of uncontrolled pressure fluctuations in 
the pump. 

But as long as the flow is governed by the linear Stokes equation (as 
it is in our case), the only effect of fluctuations of is a linear change of 
time scale. To account for this change, we plot orientational trajectories as a 
function of distance d along the trajectory of the rod (and not as a function 
of time t): 



Jo 

The instantaneous flow velocity Ux is estimated by the centre-of-mass veloc- 
ity. We assume, in other words, that the centre-of-mass of the rod is advected 




(2) 
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by the channel flow. The change of variables Q greatly simplifies the analysis 
of our results. 



3 Results 

In order to determine the sensitivity of the orientational motion, we record 
it in the following way. We first track a rod from one end of the channel to 
the other, recording the centre-of-mass position and orientation of the rod. 
Then we reverse the flow direction smoothly by reverting the applied pressure 
difference. In Stokes' equation this turns Ux to —Ux and, by Eq. ([2|, d to —d^ 
effectively reversing time. In an ideal and noise- free experiment the centre- 
of-mass and orientation are expected to retrace their trajectories. In reality, 
of course, we observe deviations which allow us to quantify the sensitivity of 
the observed dynamics to perturbations. 

Our results are summarised in Figs. |5] to [7) The data shown in the three 
figures were obtained for three different rods. Each figure contains three pan- 
els: (a) a high-resolution plot of the centre-of-mass trajectory exhibiting small 
fiuctuations Az in the z-direction, (b) the n^^-component of the orientation 
vector, before reversal (solid line), and after reversal (dashed line), and (c) 
the -component of the orientation vector. 

Consider first Figure, [s] Panel (a) shows that the fiuctuations Az of the 
centre-of-mass trajectory are very small, only on the order of 10 \im. This is 
less than one particle length, and is much smaller than the channel width 
(2500 |am). Panel (b) shows n^^, the component of orientation vector along 
the fiow direction. The particle spends most of its time aligned with the 
fiow, but regularly tumbles between Ux = I and Ux = —1. Upon reversal, 
the orientational trajectory is initially similar, but drifts out of phase after 
several fiips. The z-component of the orientation vector in panel (c) is, of 
course, close to = while the particle is aligned along the fiow. However, 
during a fiip the values of determine which particular orientational path 
the rod takes. We see that the shape of the n^-curve during fiips is similar 
for several subsequent fiips. This indicates that the orientational trajectory 
is approximately periodic. 

Fig. [6] shows a second example. It is similar to Fig. [5j but corresponds 
to a different orientational trajectory. The other main difference is that the 
period time of the tumbling motion is shorter. 

Fig. [7| shows an orientational trajectory that exhibits large deviations 
after reverting the fiow. The orientational motion is at first retraced, but 
after the rod has retraced a full fiip, deviations in the orientational motion 
become noticeable and continue to grow rapidly. 

4 Discussion 

Figs. |5] to [7| show both translational and orientational trajectories of rods 
advected in a channel fiow. The translational trajectories indicate that the 
experiment is successfully controlled: the rods stay on straight streamlines 
during the whole experiment in Figs. |5] to j?! On a finer scale we see that 
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Fig. 5: (a) Fluctuations Az of the 2:-coordinate of the centre-of-mass po- 
sition of a rod advected in the microchannel flow from left to right (solid 
black line), reversed flow (dashed black line). Here d denotes the distance 
covered in the x-direction, see Eq. ([2|. (b) n^; -component of orientation 
vector, (c) n^-component of orientation vector. Rod #1 (length 19 jam, see 
Fig. [4]). 



the streamline followed by the centre-of-mass is slightly curved, on the order 
of 10 \im. This is probably due to imperfections in the channel walls. After 
reversal, the center-of-mass trajectories of the rods follow the streamlines 
very well, with only minor fluctuations on the order of one micrometer. 

Our main results are the orient ational trajectories of the rods shown in 
Figs.[5]to[7[ As expected we observe tumbling motion, and the results include 
both periodic tumbling (Figs. 5|6) and aperiodic tumbling (Fig. [t]). Different 



realisations of the experiment result in different tumbling periods. Moreover, 
we observe a weak phase drift between trajectories before and after reversal. 

We discuss these observations by first briefly recalling the hydrodynamic 



Jeffery (1922) theory and its implications for both axisymmetric and triaxial 
rods. Small axisymmetric rods suspended in a shear flow tumble periodically. 
More precisely, the orientation vector follows one of infinitely many possible 
periodic Jeffery orbits. Which particular orbit the rod takes is determined 
by its initial orientation. For most Jeffery orbits, the rod tumbles: it stays 
aligned with the streamline for most of the time, but at times it rapidly flips 
by 180 degrees. When the initial orientation of the rod is very close to the 
vorticity direction, the tumbling is less pronounced. The corresponding type 
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Fig. 6: Same as Fig. [5] but for rod #2 (length 22 |j.m). 
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Fig. 7: Same as Fig. [s] but for rod #3 (length 30 |xm). 
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of orient at ional motion is often referred to as "kayaking" in the literature, 
owing to its similarity to the motion of a kayak paddle. 

For Jeffery orbits, the tumbling period does not depend on the initial 
orientation of the rod. The period for a rod of aspect ratio A suspended in a 
shear flow of strength s is given by 

^Jeffery =7r — — — . [6) 

As mentioned in the introduction, asymmetric particles also tumble, but not 



necessarily in a periodic fashion (Hinch and Leal, 1979 Yarin et al. , 1997). 

We begin the discussion of Figs. [5j to [7| by noting two consequences of 
Eq. p|. First, this equation shows that the different tumbling periods in 
Figs. ]5| and [6] can result from differences in aspect ratio or shear strength: 
larger aspect ratios correspond to longer periods, stronger shear to shorter 
periods. The lengths of the particles shown in Figs. |5] and |6] are approximately 
the same. But the resolution of our microscope does not allow us to precisely 
determine the thickness of the two rods. We cannot exclude that the different 
tumbling periods are in part caused by different aspect ratios of the two 
particles. 

Second, we discuss the phase drift between the forward and backward 
orientational trajectories seen in Figs, [s] and [g] in terms of Eq. ([3|. We argue 
that the phase drift is due to fluctuations in the shear that are caused by small 
fluctuations of the ^/-coordinate of the centre-of-mass position. As explained 
in Sec. |2.1[ the shear is a function of the ^/-coordinate. If we assume that the 
centre-of-mass drift in y is similar in magnitude to that in z (on the order of 
a few micrometers), we can find a lower bound for the expected phase drift 
as follows. According to Eq. ([3|, the period is inversely proportional to the 
shear strength. It follows that the relative change in shear strength equals 
the relative change in time between flips, As/s = AT/T. The flow profile 
in the ^/-direction is approximately a quadratic function, Ux = —OLy^^ which 
implies that AT/T = Ay/y. Equating y with the channel depth (400 ^im) 
allows us to estimate that the drift should at least be 1 %, in agreement with 
panels (b) and (c) in Figs, [s] and |6] 

We conclude the discussion by examining the trajectory in Fig. [7| We 
see, first, that the trajectory is clea rly aperiodic , both before and after rever- 
sal. The numerical experiments by Yarin et al.| ( |1997[ ) indicate that triaxial 



rods may exhibit aperiodic dynamics even for very small deviations from 
axisymmetry (that could not be resolved by our microscope). Secondly, the 
forward and backward orientational trajectories in panel (c) of Fig. [t] sepa- 
rate rapidly. In fact, the trajectory appears at first to reverse perfectly but 
then the difference between forward and backward orientational trajectories 
grows substantially. These observations suggest that Fig. [7| corresponds to 
chaotic tumbling of a triaxial particle. 



5 Conclusions 



We have designed a microfiuidic setup with video microscopy and tracking 
software to measure the translational and orientational trajectories of micro- 
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rods advected in microchannel flows. The experiments presented here demon- 
strate the level of control and accuracy of the current setup. We observe both 
periodic and aperiodic (and possibly chaotic) orientational dynamics. 

Further work on the experiment aims to improve the efficiency and capa- 
bility of the current setup. For example, we plan to install an optical tweezer 
in order to control initial orientations of the rods. This would make it possible 
to observe different orientational behaviours (such as those shown in Figs, [s] 
and [t]) for the same rod. 
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